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Abstract 

We discuss the case of a Markovian master equation for an open system, 
as it is frequently found from environmental decoherence. We prove two 
theorems for the evolution of the quantum state. The first one states 
that for a generic initial state the corresponding Wigner function be- 
comes strictly positive after a finite time has elapsed. The second one 
states that also the P-function becomes exactly positive after a deco- 
herence time of the same order. Therefore the density matrix becomes 
exactly decomposable into a mixture of Gaussian pointer states. 

PACS numbers: 03.65Bz, 05.60.Gg 



1 Introduction 



The study of Markovian open systems is of general interest. One of the rea- 
sons is that the coupling of open systems to their ubiquitous environment often 
leads to master equations which are local in time |T[. The interaction of dust 
particles with air molecules or radiation, for example, delocalises any inter- 
ference terms into correlations with the environmental degrees of freedom on 
a short decoherence timescale. Thereafter, the dust particle can be perfectly 
described by a Markovian master equation for its density matrix p(t): 

Such an equation results frequently from an interaction with the environment 
in cases where friction is negligible • The strength of the coupling is given 
by the parameter D. The first term in (|l|) would lead to unitary spreading, 
whereas the second term would lead to nonunitary localisation. For a wave 
function of characteristic width a, these effects act on timescales ma 2 and 
1/Da 2 , respectively. Both effects are thus balanced for an "equilibrium width", 
given approximately by 

a = (Dm)- 1 /' , (2) 
see |2|, [l|]. The corresponding timescale is 



t = ^m/D (3) 

and will set the timescale for decoherence. 

The oldest way to elucidate a quantum state in terms of a pseudo-classical 
distribution of phase-space variables T = (x,p) is due to Wigner. Instead of 
the density matrix, one discusses the Wigner function 

W(T) = f(x - \r\p\x + ±r)e ipr dr (4) 

in order to study aspects of an open system . Normalization holds with the 
choice dT = dpdx/2ir of the phase-space volume element. The Wigner function 
W(T; t) corresponding to ([!]) satisfies the Fokker-Planck equation 

dW p dW Dd 2 W 
dt m dx 2 dp 2 

In a general situation, the Wigner function is negative in some regions of phase 
space. For this reason it cannot be regarded as a probability distribution. If 
decoherence occurs, the general expectation is that these negative parts are 
smoothed out in the course of time. Many examples support this expectation. 
In Sect. 2 of this paper we shall prove a much stronger statement: After a 
certain decoherence time, the Wigner function becomes strictly positive. This 
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is very different from the behaviour of the density matrix whose nondiagonal 
terms (describing interferences) become zero only asymptotically, i.e., they 
remain nonzero at any finite time. 

Through the process flip or, equivalently, fl^) the position basis is distin- 
guished as the preferred basis with respect to which no interferences can be 
observed and which remains a robust basis in time. Such a classical basis 
is usually called the pointer basis. It is the basis which remains most sta- 
ble against the influence of an environment. How can the pointer basis be 
determined? We have shown in Q that this can be achieved by three dif- 
ferent methods which all lead to the same results. The first method invokes 
the principle of "Hilbert-Schmidt robustness" stating that the pointer states 
mimic the local nonunitary evolution as closely as possible with respect to the 
Hilbert-Schmidt norm. A unique set of Gaussian pointer states has thereby 
been found. The second method demands that the local production of entropy 
be minimal ("predictability sieve"). Again, this has led to a set of unique 
Gaussian pointer states with practically the same width as the previous ones. 
The last method invokes quantum state diffusion and leads again to Gaussian 
pointer states with the same width. 

In Sect. 3 we introduce the overcomplete set of Gaussian pointer states 
Then we define additional phase-space distributions: The (generalized) Q- and 
P-functions are related to the Wigner function by Gaussian coarse-grainings. 
In the discussion of the local entropy production || we have already used 
the important fact that the density matrix can be decomposed exactly into a 
mixture of Gaussian states after a finite decoherence time. In Sect. 4 of the 
present paper we shall prove this theorem - the positivity of the P-function 
- which, in fact, holds for arbitrary initial states. The technical details are 
relegated to Appendices A, B, and C. 

2 Strict positivity of the Wigner function 

We shall now consider the Wigner function, W(T;t), of our open quantum 
system, obeying Eq. (|5|). We shall prove the theorem that for any initial state, 
W will become positive after a finite time t^, i.e., 

W{T;t)>0, t>t D . (6) 

W can of course be positive even earlier; the theorem states that it cannot 
be negative later. The Wigner function corresponding for a Gaussian wave 
function, for example, stays positive for all times. In fact, as shown e.g. in [|J, 
Gaussian states are the only pure states that lead to a positive Wigner function. 
It must, however, be emphasised that positivity is only a necessary requirement 
for classicality, not a sufficient one. Squeezed Gaussian states of harmonic 
oscillators, for example, are genuine nonclassical states, but correspond to a 
positive Wigner function. 
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Since Fokker-Planck equations like (^) preserve the positivity of distribu- 
tion functions, one might guess that the set of positive solutions attracts the 
indefinite ones, so that any indefinite distribution becomes positive in the 
course of time. This is, however, not true. Eq. (||) is a linear equation. Since 
it also allows solutions that remain negative, the superposition principle will 
allow indefinite solutions at all times. 

The proof, therefore, has to rely on specific properties of the Wigner func- 
tion. It is know that W possesses, in fact, very special features, see e.g. P. 
For example, it obeys \W\ < 1/tt. Moreover, its negative parts are always 
restricted to small regions in phase space. 

To be concrete, in our proof we shall make use of two properties for the 
Wigner function. The first one exploits its connection with the Q-function (see 
e.g. ||). Introducing the normalized Gaussian 



g(T;C) = ^exp 



x 2 + p 2 
2C 



(7) 



its convolution with the Wigner function, g * W, is always positive provided 
C > 1/2 (and may be indefinite otherwise). In the marginal case C = 1/2 this 
convolution gives the Q-function which yields the probability distribution for 
finding the coherent states in the density operator p, and is therefore manifestly 
positive || . This is why the negative regions for W are so restricted: the con- 
volution of W with the Gaussian g(T; 1/2) would not yield a positive function 
if W contained negative regions that are spread out over regions with areas 
much bigger than about 1/2. Since the case C > 1/2 represents a stronger 
coarse-graining of W than C = 1/2, it is clear that this positivity remains 
true. 

The second property of the Wigner function, which will serve as an ingre- 
dient in our proof, is its invariance under linear canonical transformations. To 
be precise, if we make such a transformation for both the quantum operators 
and the classical canonical variables, one has 

(x,p) -> {x',p') => W(x,p) -> W'(x',p') = W(x,p) . (8) 

This feature is not obvious from the usual expressions defining W. There 
exists, however, an alternative equivalent form given by ||: 

W(x,p) = ti[p{5(x - x)5(p-p)} sym ] , (9) 

where we refer to a simple symmetrisation process regarding the order of x 
and p: 

{f6} sym = -(td + Of), for f = (x,p) and for all 6. (10) 

This ordering is explicitely invariant for linear canonical transformations. The 
invariance of the Wigner function follows immediately. Note, however, that 
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the coarse-graining is not invariant and coarse-grained Wigner functions, Q- 
functions in particular, will be non-invariant even for linear canonical trans- 
formations. In the simple case of the transformation x' = ax,p' = a~ 1 p, and 
assuming a pure state with wave function ip(x), the invariance can be seen 
directly from the standard integral expression. Using 

W) = ^ (-) , (ii) 



one has 



W'(x',p') = - f dy'exp{2ip'y')4j*'{x' + y'U'(x' - y') , (12) 

71 J 



where W denotes the Wigner function with respect to the same state in the 
transformed basis, and therefore 

W'(x',p') = W(x,p) . (13) 

Using this invariance, one can extend the previous coarse graining (|7|) to gen- 
eral Gaussians with correlation matrix C: 



ff (r;C) = |C|-2exp 



rp 1 ' 

-r T — r 

2C 



(14) 



(When vectors, V stands for column and T T for row vectors, respectively.) 
Application of a linear canonical transformation to the convolution g * W 
rendering C = J |C|I then demonstrates that the sufficient and necessary 
condition for the positivity of the coarse-grained Wigner function reads 

g{T-C)*W(T) >0 iff |C|>l/4. (15) 

With this lemma, the proof of the theorem @ becomes straightforward. 
The Fokker-Planck equation (^) imposes a progressive Gaussian coarse-graining 
([14]) on the initial Wigner function: 

W{T- 1) = g{T; C w (t)) * W{x - pt/m,p; 0) , (16) 

where the time-dependent correlation matrix of the coarse-graining is 

'iff*?")- m 

as can be found from Eqs. (^|, ^4| . ^6|) . The determinant yields 

ICVWI = ^ . (18) 
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It follows from the condition QT5| ) that the Wigner function is indeed positive 
for 

— > 3 1/4 » 1.32 , (19) 

which is of the same order as the decoherence timescale of Sect. 1. This 
completes our proof. 

We would like to mention that, based on more complicated mathematics 
0, the lemma (|l^) has been stated earlier ||. Similarly to the spirit of [10] 
and the present work, it was also suggested in |7|, || that particular quantum 
features of open systems would disappear after an estimated finite time. 



3 Pointer states, Q- and P-functions 

We are going to consider the overcomplete set of normalized pure Gaussian 
pointer states |T) obeying 

/ |r><r|dr = / . (20) 

In Appendix A we discuss the position representation of these states. The 
overlap of two pointer states is a Gaussian: 



l<r|r'>| 2 = exp 



ir _ r y (r _ r) 

41^1/4 



(21) 



The matrix C1/4 is positive and has determinant 1/4. This value for the 
determinant makes the above Gaussian normalized and assures the consistency 
of the completeness relation ( f20| ) with the normalization of the pointer states. 
In fact, C1/4 can be directly identified as the matrix of quantum uncertainties 



of the canonical pair of operators in the Gaussian pointer states, cf. in 

Appendix A. We can calculate C1/4 for simplicity in the fiducial state |T = 
(0,0)) = |0>: 

Ci/4 = (0|ff T + h.c.|0) . (22) 

Given the above overcomplete set of Gaussian states, one can introduce the 
generalized Q- and P-functions, related to the density operator respectively by 

Q(T) = (T\p\T) (23) 

and 

p = Jp(T)\T)(T\dT . (24) 

The Q-f unction is the probability distribution of the value T. In a generalized 
quantum measurement it can be inferred from the positive operator valued 
measure (see e.g. |§) formed by |r)(r|dr. The P-function has a different 
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meaning. When we expand the density matrix as a sum of the pointer states, 
the weighting function is called the P-function. While the Q-function is non- 
negative by construction, the P-function may be indefinite (even ill-defined) 
for generic states. The Q- and P-functions are related to the Wigner function 
by the same Gaussian coarse-graining but in opposite senses 0: 

w/(r) = <7(r ; c 1/4 )*p(r) , (25) 

Q(T) = g(T; C 1/4 ) * W{T) . (26) 

The correlation matrix C1/4 of coarse-grainings is the one that appeared ear- 
lier as the matrix of canonical quantum uncertainties in the Gaussian pointer 
states, see (p2|). All sets of Gaussian pure pointer states are classified in Ap- 
pendix A. Details of derivations for (|25|) , (f2~6|) are given in Appendix B. 

The Q-function satisfies a Fokker-Planck equation which is, according to 
(p6j), the coarse-grained version of the Fokker-Planck equation (|J) for the 
Wigner function. For the evolution equation of the P-function we are going to 
present a more direct derivation in Appendix C. 



4 Strict positivity of the P-function 

In Sect. 3 we have expanded the density operator as a weighted sum of Gaussian 
pointer states, see fl24|) . The weight function is called P-function and it is 
indefinite for a generic quantum state p(0). We shall now prove that, due to 
the open system dynamics ([!]), the P-function becomes exactly positive, 

P{T;t)>0, t>t' D . (27) 

The density operator p(t) can thus be decomposed exactly into a statistical 
mixture of Gaussian pointer states |T) (see ( P2"| ) below) after a finite decoher- 
ence time t' D has elapsed. We already emphasise here that this holds for a 
generic p(0) (not necessarily Gaussian). This theorem generalises the corre- 
sponding statement made in [pl| for a single choice of Gaussian pointer states 



(i.e. of C1/4) as well as the asymptotic statements proved in [IT], [12 

The proof will be reduced to the lemma (|T5|) used in the proof of the 
positivity of the Wigner function, see (§) in Sect. 2. Consider the solution ( |16"D 



for the Wigner function and substitute the expression ( p5|) into its left-hand 
side. As for the right-hand side, assume that enough time has elapsed so that 
Cw{t) — C1/4 is a non-negative matrix. Then the convolution factorizes as 

g(T; C w (t)) = g(T; C 1/4 ) * g(T; C w (t) - C 1/4 ) (28) 

and can be substituted into (|T6"p. We obtain 

g(T; C 1/4 )*P(r; t) = g(T; C 1/A )*g(T; C w (t)-C 1/4 )*W(x-pt/m,p; 0) . (29) 
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The identical convolutions on both sides cancel each other and leave us with 
the explicit solution for the P-function as the coarse-grained Wigner function. 
The lemma fll5|) tells us that a generic initial P-function becomes non-negative 
after a time t, 

P(r;i)>0 iff |CV(t)-C 1/4 |>l/4. (30) 

Calculating the determinant from ( |T7| ) and fl3T)| ) below, this condition leads to 
a cubic equation for t, yielding the numeric estimate 

- > 1.97 . (31) 
to 

As expected, therefore, the decoherence timescale t' D coincides approximately 
with to- This completes the proof. 

One might wonder how much of our results do depend on our particular 
Markovian master equation ([I]). We expect that our proofs of positivity can 
be extended to all master equations of this kind which are at most quadratic 
in position and momentum, e.g. for the damped harmonic oscillator. As far as 
corrections to the Markovian approximation are concerned, we do not expect 
that they render our results - exact positivity of Wigner and P-function after 
a finite time - obsolete; because of their exponential smallness they should not 
spoil an exact positivity. Our results should thus be stable with respect to 
post-Markovian corrections, but a formal proof would go beyond the scope of 
this paper. We hope to return to these issues elsewhere. 



A Classes of Gaussian pointer states 

We choose the following form for the wave functions {q\T) of pure Gaussian 
pointer states 

MQ) = (a R /2ir) 1/4 ex P (-a(q-x) 2 /4 + ip(q-x)) . (32) 

This is the most general form of a normalized Gaussian wave packet shifted 
and boosted uniformly from the fiducial state for ipo(q) for V = (0, 0), 

MQ) = (« R /2vr) 1/4 exp(-«g 2 /4) . (33) 



The correlation matrix C1/4, defined in (22|), takes the following form: 

(01 ( * (xp + h.c.)/2\ J_( 1 -otj/2 \ 

{[)l {(xp + h.c.)/2 f ) |U) " aR { -on/2 \a\ 2 /A)- {6i) 

Indeed, its determinant is 1/4 for all a. In [||] we have found a distinguished 
value of the complex parameter a for the class of robust pointer states, given 
by 

a = an + iaj = (1 — \)\/2Dm . (35) 
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For this value ( pip gives 




This is being used in the proof for the positivity of the P-function, Eq. (j30|). 
The proof can be performed for any a (it always leads to a cubic equation), 
but the basic object is the matrix Cj/4 with determinant 1/4. 

It is worthwhile to add the following. Let us assume that we had con- 
structed no pointer states first. Instead, assume that we had taken an arbi- 
trary real positive matrix of determinant 1/4 to perform a coarse-graining on 
the Wigner function like in (p6|). This matrix would then define a complex pa- 
rameter a as on the right-hand side of (0) , and we would be able to construct 
a unique overcomplete set of Gaussian wave functions fl32|) whose Ci/4-matrix 
is just our chosen one. Thus it would turn out that our coarse-grained Wigner 
function were just the corresponding Q-function which is always positive. In 
this way we have obtained an alternative proof of the theorem (0) of Sect. 2. 



B Symplectic Fourier transform 



This appendix is devoted to the proof of (|25|) and (|26|) . To facilitate the 
mathematical derivations we use the Fourier representation. For instance, the 
Fourier transformed Wigner function reads 



W(T) = / W(T) exp 



-if T r 



dr 



(37) 



where f = (p, —x). We use the symplectic product T T r = px — xp of the orig- 
inal and the Fourier-transformed variables in fl37f) . In Fourier representation, 



the coarse-graining relations (p5| , |26|) reduce to algebraic ones: 



W(T) = exp 



Q(T) = exp 



2 

1~ 



r T c 1/4 f 



r J c 1/4 r 



p(r) 



W(T) . 



(38) 



(39) 



We are going to prove the first relationship (the second one can easily be 
proven along the same lines). From (0J) and (^) we obtain 

W(t) = J e ii,x (x + ±x\p\x - \x)dx . (40) 

Substituting the P-function expansion (|24[) of p on the right-hand side yields 

W(t) = f e~ ipx (x + i5|r')(r'|a; - ix)P(r / )dr' . (41) 



s 



Invoking the expression (|32[) of the pointer state wave functions we can perform 
the above integral, yielding 



W(f) = exp 



a R x 2 



2a 



R 



\p-\afx) 2 P(T) 



(42) 



The quadratic form in the above exponent can be written as — ±r T Ci/4r. This 
completes the proof. 



C 'Fokker- Planck' equation for the P-function 

In this part we return to the notations of our recent work j3| and derive directly 
the evolution equation for the P-function. We introduce the projectors on the 
states i[> r (q) = (q\T), P(T) = \T)(T\. We also denote the P-function P(T;t) of 
the open system by f(T;t). Application of (|1|) to (f24j) yields 

= / /(T; t)CP(T)dT = J f(T; t)P(T)dT . (43) 

To derive from this the equation for /, we first observe that the action of C 
on P(T) in coordinate representation yields 

{q\CP\q) = — < — i«r — ip \a*(q' — x) + a(q — x)] 
2m [ 

^)\q'-x) 2 +(f)\q-x?}(q\P\q') 
-§(q-q') 2 (q\P\q') ■ (44) 

Using fl32|) in the identity (q\P(T)\q') = V'r(<?)V'f W)i ^ ne above action can 
be expressed as an operator containing derivative terms in x and p, acting on 
P. After partial integration in (|43|) one then finds the following equation for 

/, 

= -^d x f(T; t) + ~ [D pp d 2 pp + D xx dl x + 2D px d 2 px ] f(T; t) , (45) 

where the elements of the diffusion matrix are given by 

n _ ( D xx D xp \ _ ( -aij man \a\ 2 /Ama R \ , , 

U ~ \ D px D pp )-{ \a\ 2 /Ama R D ) ' ^ 

Eq. ( ^§1) can be interpreted as a Fokker-Planck equation for f(T,t) provided 
the diffusion matrix D is non-negative. As usual, the first term in fllS]) will then 
describe a drift according to the free-particle dynamics, while the second term 



9 



will describe a diffusion of the state of the system over the pointer states P(T). 
In we have implemented the principle of minimal local entropy production 
by minimising the width of the Gaussian pointer states. This leads to a value 
for a of the order of (|35|). We have shown that the relation C = Dig/2 holds 
for the distinguished value fl35p of a. 

To find a formal solution of (^Sj), we use the Fourier representation f(T; i). 
Eq. (|45|) then leads to the following evolution equation for the Fourier compo- 
nents: 



at m 2 



r^Dr 



Km . 



The solution assumes the form 



/(f; i) = exp 



-^FD(i)f 



f(x — pt/m,p; 0) . 



(47) 



(4£ 



By substitution into (|47|) one obtains explicitly the matrix of time-dependent 
coefficients 



D(i) 



Dxx + D xp t/m + D pp t 2 /3m 2 D xp + D pp t/2 



m 



D xp + D pp t 1 2m 



D 



(49) 



pp 



This matrix becomes more and more positive for t > provided D = D(0) 
was chosen positive. One obtains the solution of the evolution equation (|45l) 
by the inverse Fourier-transform of the expression (48). It takes the form of 
the convolution 



/(r; t)=g {V; iD(i)) * f{x - pt/m, p; 0) . (50) 

The solution therefore emerges as the progressive time-dependent Gaussian 
coarse-graining of the free kinematic evolution. Note the close similarity of 
this equation with the time evolution (O) for the Wigner function. 
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